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Introduction. The paper is devoted to the study on free oscillations of the sensing element of a micromechanical R-R- 
type gyroscope of frame construction developed by the Kuznetsov Research Institute of Applied Mechanics, taking into 
account the nonorthogonality of the torsion axes. The influence of the instrumental manufacturing error on the accuracy 
of a gyroscope on a movable base in the case of free oscillations is studied. The work objective was to improve the 
device accuracy through developing a mathematical model of an R-R type micromechanical gyroscope, taking into 
account the nonorthogonality of the torsion axes, and to study the influence of this error on the device accuracy. The 
urgency of the problem of increasing the accuracy of micromechanical gyroscopes is associated with improving the 
accuracy of inertial navigation systems based on micromechanical sensors. 

Materials and Methods. A new mathematical model that describes the gyroscope dynamics, taking into account the 
instrumental error of manufacturing the device, and a formula for estimating the error of a gyroscope, are proposed. The 
dependences of the state variables obtained from the results of modeling and on the basis of the experiment are 
presented. Methods of theoretical mechanics and asymptotic methods, including the Lagrange formalism and the 
Krylov-Bogolyubov averaging method, were used in the research. 

Results. A new mathematical model of the gyroscope dynamics, taking into account the nonorthogonality of the torsion 
axes, is developed. The solution to the equations of small oscillations of the gyroscope sensing element and the estimate 
of the precession angle for the case of a movable base are obtained. A comparative analysis of the developed model and 
the experimental data obtained in the case of free oscillations of the gyroscope sensing element with a fixed base is 
carried out. The analysis has confirmed the adequacy of the constructed mathematical model. Analytical expressions are 
formed. They demonstrate the fact that the nonorthogonality of the torsion axes causes a cross-influence of the 
amplitudes of the primary vibrations on the amplitudes of the secondary vibrations of the sensing element, and the 
appearance of an additional error in the angular velocity readings when the gyroscope is operating in free mode. 
Discussion and Conclusions. The results obtained can be used to improve the device accuracy using the algorithm for 


analytical compensation of the gyroscope error and the method for identifying the mathematical model parameters. 
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Introduction. The development of high-precision micromechanical inertial sensors, including 
micromechanical gyroscopes (MMG), used to solve navigation problems and control the movement of aircraft and 


mobile robots, is an urgent task of instrument engineering [1]. The advantages of MMG include small weight and size, 
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as well as low cost compared to gyroscopes based on other physical principles. However, the main disadvantages of 
MMG are the variability of its metrological characteristics and the low accuracy of measuring the parameters of the 
angular motion of the object (angular rate and angle of rotation). The principle of operation of vibrating gyroscopes is 
based on the property of the Foucault pendulum to keep the plane of small vibrations motionless in inertial space [2]. 

Fundamentals of the theory of gyroscopes of the generalized Foucault pendulum class, which include MMG, 
are given in [2—5]. They describe various design schemes for the construction of MMG, and investigate the influence of 
instrumental manufacturing errors and changing operating conditions on the dynamics of the gyroscope. The principal 
feature of the gyroscopes of the generalized Foucault pendulum class is nonlinearity due to the finite vibrations of the 
sensing elements (SE) or physical nonlinearity associated with the features of the vibration control system [2-6]. 

Studies on the MMG dynamics and design were also published in the works of foreign authors [6—9]. For 
example, in publications [6, 9], a formula for estimating gyro drifts was obtained, based on the use of a developed 
mathematical model of motion that describes a slow change in the toroidal coordinates of the SE vibrations. In [7, 8], 
the issues on manufacturing MMG are discussed, and the equations of its small vibrations are analyzed. In paper [7], the 
equations of the MMG motion with angular (R-R-type) and linear (L-L-type) oscillatory types of the SE movement are 
compiled. In the above paper, a comparative analysis of the dynamics of such devices is carried out within the 
framework of linear models, and recommendations are given on the selection of MMG parameters, based on the 
conditions for increasing sensitivity and ensuring the required bandwidth, as well as the requirements for the linearity of 
the scale factor. 

When designing MMG, the developers tend to use the phenomenon of internal resonance in the system, due to 
the combination of the natural frequencies of the SE vibrations [3, 4]. However, it is noted in [7, 8] that errors in the 
manufacturing technology, unknown and unpredictable deviations of structural elements from the design positions 
cause additional errors in the device measurements. 

To improve the accuracy of the measurement of the angular rate of the MMG, the objective is set: to study free 
oscillations (in the absence of control) of the R-R-type MMG SE, taking into account the effects arising from the 
nonorthogonality of the torsion axes. This defect appears due to the imperfection of the manufacturing technology of 
the device. The tasks of developing a new mathematical model of the MMG dynamics considering the nonortho gonality 
of the torsion axes, evaluating the device drift, and describing the effect of the nonorthogonality of the torsion axes on 
the dynamics of the MMG SE, are set. 

Materials and Methods. A model design of an R-R-type vibrating MMG — a design with an intermediate 
frame in accordance with the classification from the source is considered [3]. The kinematic scheme of the gyroscope 


(Fig. 1) is implemented in the form of a two-degree gimbal of the SE. 
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Fig. 1. Design diagram of the device: 1) base (body); 2) intermediate (external) frame; 3) sensing element consisting of balanced 
plate and inertial mass; 4) torsion bars 


To describe the SE position, we introduce coordinate systems (Fig. 2) associated with: the device body — 





232 OXYZ ; with the external frame of the elastic suspension of the gyroscope — Ox,y,z,; with the balanced plate 
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Oxyz . Moreover, OZ is the axis of sensitivity of the gyroscope, and the coordinate system Ox, y,z, differs from the 
system Ox,y,z, by turning by a constant angle of nonorthogonality of all torsions. In the presented systems, the origin 


of coordinates corresponds to point O and is located in the geometric center of the balanced plate. 





Fig. 2. Coordinate systems 


In the system Ox,y,z, we will set the axial moments of inertia of the intermediate frame Jj,, J,, Jj,, and in the 
system Oxyz — the axial moments of inertia of the SE Jp,, J2,, J2,. Note that in this paper, the axes of the coordinate 
systems Ox,y,z, and Oxyz are considered the main central axes of inertia of the intermediate frame and SE, 
respectively. 

When modeling the motion of the SE, the assumption is made that the torsion design provides infinite bending 


stiffness. The position of the SE relative to the base of the MMG is described by two generalized coordinates — angles 


a and £, as well as a small constant angle 6 , that characterizes the nonorthogonality of the torsion axes (Fig. 2). The 
relative position of the coordinate systems is determined by a sequence of elementary rotations: 


OXYZ —*-> Ox, y,z, > Ox, yz, > Oxyz , 
x 2) 2 


where under each arrow, the axis is indicated around which there is a counterclockwise rotation by the angle indicated 
above the corresponding arrow. 

We will set up the equations of the dynamics of the MMG SE in the form of Lagrange equations of the 2nd 
kind [10, 11]: 
ae, eee " 
dt\0a) 0a da adt\op) BB ap’ 
where L = 7—JI — the Lagrange function; 7 and /7 — kinetic and potential energies of the system, respectively; ® — 


dissipative function that characterizes the loss on internal friction. The expressions for these values have the form: 


2x "x 2y"y 22°" 2 


2 
on14ctstag%, naleatades ° 
7 ye i 


where d,, d oo friction coefficients; c_, Cz — torsion stiffness coefficients. 


The expressions for the projections @,, ®,, ®, of the SE angular rate on the movable axes x, y, z have the form: 
o, =acosBcosd —-QcosasinB+QsinacosPsino, 
@, =B+Qsinacosd-asino, (3) 


o, =asinBcosd+QcosacosB+QsinasinBsino. 
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Given the smallness of angles a , B and o, the trigonometric functions in expressions (2) and (3) from these 
angles can be replaced by Taylor series expansions, limiting the summands to the first order of smallness. Then we 
obtain the equations of small vibrations from the equations of motion (1), written with account for the expressions (2) 
and (3). 

Using the Lagrange formalism [10] in the case of a constant angular rate of the base, we obtain the equations of 
small oscillations of the SE, written up to the terms of the first order of smallness in the form: 

é+ota=j06-°a+5f, +028 =-j,04-““ f+ 04, (4) 
QO, jh QO, 


where the following notation is introduced (similar to how it is performed in paper [12]): 








- J. +J),, —J), . J 5. +J,, —J), Cy 
J ~~ J J 2 J. = J 9 QO, a 9 
1x as 2x 2y Ji + Jy. 
®, = Cp om = d, QO." = dp 
B = p) a. p) B ~ 
2¥ 0, (Ji, +J,,) OJ 5, 
Here, j,, j, — dimensionless moments of inertia of the elastic suspension; @,, Og and Q,, Q, — natural 


frequencies of vibrations and Q-values at corners @ , B , respectively. 

When deriving the oscillation equations (4), the angular rate of the gyroscope body © was considered small 
relative to the natural frequency @,, T.e. IQ <<q@,, and angle 6 was also assumed to be a small value, ie., 06 <<1. 
Note that 1n equations (4), the terms are dropped due to the presence of geometric nonlinearity of the MMG. The 
influence of the nonlinearity of the geometry of the SE motion on the dynamics of the R-R-type MMG 1s described in 
monograph [5]. 

Taking into account that the right-hand sides of equations (4) are small perturbations, i.e., @+@,a@ = O(e) , up 


to the terms of the first order of smallness, we can write: @ = —0. a + O(e) . Thus, the second derivatives of angles @ 


and $6 are excluded from the right part of equations (4). 
We consider the case of an isotropic elastic suspension, 1.e., the equality of natural oscillation frequencies and 
equal Q-factors: 
O, = Wg =, 0, =9, =9, 


where @, — the characteristic value of the natural frequency of vibrations; Q — the characteristic value of Q-factor. 


It should be noted that the case of a difference in quality (O, # O, ) and a small difference in frequency 


(a, # Op | under studying free oscillations of the MMG SE is considered in paper [12]. With the entered notation and 


the accepted assumptions, we write down the equations of the SE motion with an accuracy of the first-order terms of 
smallness in a dimensionless form: 
&+o.a = j, OB-O'o,a- A Sop, 
J (5) 
B+op =—j,Q24-O'w,B-Sa.a. 

Note that the system of equations (5) is reduced to the standard form of writing a regularly perturbed system of 
differential equations with one fast angular variable [13, 14]. One of the most common ways to find solutions to 
regularly perturbed systems is to use asymptotic motion separation methods [13-16]. 

The solution to nonlinear equations (5) is obtained using the Krylov-Bogolyubov averaging method [14], and 
we will use Van der Pol variables as slowly changing variables [13] p,, 41, D2, g2: 

A = p,Sin(@t)+¢,cos(@t),  G&=,p,cos(@t)—ap gq, sin(@,f), 


B= D> sin (@)t) + qd cos(@,t), p = Wo P2 COS(@t)— Oy, sin (@)f). 
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Using the averaging procedure [15, 16] over an explicitly incoming time, we obtain an averaged system of 


differential equations solved with respect to the derivatives of slow Van der Pol variables: 














1 j,Q j,6 1. j,0 
pi =--O'p,+2—p,-+-¢,, gq =--0'¢,+ + p,, 
2 2Q) 2). 2 27, (6) 
1 4 Jo? O ; 1, 2 O 
S = — = = i = = +-— F 
P> 5¥ D> Wi, D 54 q) 5¥ q) 2M, 5 Pr 


The dash in equations (6) indicates the differentiation in dimensionless time tT = @t . 


The resulting model in the form of linear differential equations describes the free oscillations of a gyroscope 


SE on a movable base. The solution to the system of equations (6) can be written as: 


i 3 cos (fv? +777) + - we at car inl 


(7) 


= exp 


+V 
)=ex(-5] Dio cos (v? +7 74 A (7Pa0 +¥40) et sin( (v7 +7 malt 
20 Jr lv? yy? 
lis i (vp. + 
)=e0{ - 5 | 5] Ds cos( (Vv? +7 7 | Ja (YP +740) ee sin (v? +y °z) |, 
20 j \v 49° 
i (v 
)=er0{ - 5 Gi cos| (WV? +y 7) hr (V4 = 720) ee sin lv? +7 °r| 
20 ji |v’ +? 
where pi. = P, (0), Pry = Pr (0), do = (0), Y9 = 42 (0) — initial conditions; dimensionless angular rate of the device 


= exp 


base v=, | ee: i ( 2@, ) ; parameter characterizing the nonorthogonality of the torsion axes yv = Liv O i (2 qs , 


The second terms in formulas (7) characterize the cross-influence of primary vibrations on secondary 
vibrations and vice versa. Note that in the case of orthogonal torsion axes, when /” = () , solution (7) coincides with the 
results of paper [5]. 

The obtained analytical solutions (7) of the oscillation equations are of interest for the development of methods 
for identifying parameters, as well as predicting the gyro drift and considering it when using the method of algorithmic 
error compensation. 

Research Results. To validate the developed model, we compare the simulation results calculated from 
formulas (7) and experimental data. The measurement information was obtained using an observation system. As the 
measurement information of electrostatic sensors, we have Van der Pol variables p,, q), po, go. 


In the experiment, a sample device with the following parameters of a mathematical model with a fixed base 
(v=0): O=3856, j, =j,=1, y=0.2-10~,as initial conditions for the Van der Pol variables, values equal to the 
measurements at the initial moment of time were selected: 
1) =13.467-10°, 9, =20.429-10°, p,, =0.787-10°, gq) =1.172-10™. 
Parameter y value corresponds to the angle of nonorthogonality of torsion axes 6 equal to one angular 
second. A graphical representation of the dependencies of Van der Pol variables p,(z), 9,(7), p.(7), 9.(7) that 


slowly change over a dimensionless time 1s shown in Fig. 3. 
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Fig. 3. Dependences of Van der Pol variables: 1 — simulation results; 2 — experimental data 


The graphs (Fig. 3) show a significant coincidence of the dependencies for variables p,, g,, obtained from 
the simulation results with the experimental data. The dependences for variables p,, gq,, obtained from the simulation 


results, are in qualitative agreement with the experimental data, and the observed small quantitative deviations may be 
due to nonlinear effects, such as the nonlinearity of the geometry of the SE motion [5], or the phenomena of variability 
of quality, variability of frequency, and the error of the inertial mass displacement [12]. Consideration of nonlinear 
effects affecting the dynamics of MMG in the construction of mathematical models of SE oscillations increases the 
accuracy of micromechanical sensors as part of the inertial navigation systems [2]. 

The gyro drift due to nonlinear effects and other instrumental errors will be estimated using the auxiliary 
functional J [5, 6, 9]: 


_ 24 IIs (4% 7 DiDs) 
ary a ay er 8) 
Js (4, T Py )-A (q; T Py ) 
which is related to angle @ via the relation: 
1 
0= 5 aretan(/) . 
Moreover, this parameter is proportional to the integral of the angular rate: 


0= MAP fa(n ar, 
0 


Using formula (8), taking into account solution (7), it is possible to estimate the gyroscope drift associated 
with the nonorthogonality of the torsion axes, which arose due to the imperfection of the manufacturing technology. 
Figure 4 shows the dependences of functional J on the dimensionless time according to the results of the experiment 
and the calculation according to formula (8). 
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Fig. 4. Functional dependence / (z) : 1 — simulation results; 2 — experimental data 


The proximity of the drift estimate (Fig. 4) and the Van der Pol variables (Fig. 3) obtained analytically with 
experimental data characterize good accuracy of the constructed model, especially if we take into account the fact that it 
neglected nonlinear effects, as well as the phenomena of variability of frequency and variability of quality. Despite 
these neglections, the model provides building methods for identifying parameters, with the help of which it is possible 
to clarify the dependences obtained under modeling. The application of methods for identifying the parameters of a 
mathematical model will cause an increase in the accuracy of the MMG in the forced oscillation mode, which is the 
operating mode of gyroscopes. 

Discussion and Conclusions. A new mathematical model of the R-R-type MMG for the mode of free 
oscillations of the SE is constructed. The model takes into account the nonorthogonality of the torsion axes, which 
arises as a result of the technological impossibility to ensure high accuracy of the device manufacturing. The formula 
for estimating the precession angle with a movable base of the device is obtained. Through comparing the simulation 
results and the experimental data, the validation of the MMG mathematical model was carried out. It is shown that the 
nonorthogonality of the torsion axes causes a cross-influence of primary vibrations on the magnitude of secondary 
vibrations and vice versa. The research results can be used in the algorithm of analytical compensation of the gyroscope 
error to increase the accuracy of the MMG. 
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